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Abstract 

For a non-self-interacting Bose gas with a fixed, large number of particles 

confined to a trap, as the ground state occupation becomes macroscopic, the 

condensate number fluctuations remain microscopic. However, this is the 

only significant aspect in which the grand canonical description differs from 

canonical or microcanonical in the thermodynamic limit. General arguments 

and estimates including some of the vanishingly small quantities are compared 

to explicit, fixed-number calculations for 10 2 to 10 6 particles. 
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I. INTRODUCTION 



Large fluctuations are a salient feature of the thermal behavior of systems of bosons. 
For example, if n is the mean number of non-interacting particles occupying a particular 
one-particle state, then the mean-square occupation fluctuation is n{n + 1). This is easily 
derived in the grand canonical picture by considering diffusive equilibrium with a particle 
reservoir characterized by a chemical potential JTJ. If, however, the system has a fixed total 
number of particles, N, confined in space by a trapping potential or container, then at low 
enough temperature T or fixed total energy E when a significant fraction of N are in the 
ground state, such large fluctuations are impossible. No matter how large N, this aspect of 
the grand canonical description cannot be even approximately true. This paper addresses 
what does happen for fixed large N as iV —>■ oo. 

A decades-old answer to this question is that any interaction between the particles would 
eliminate such large fluctuations, even in the presence of a chemical potential. With a 
weak inter-particle interaction and a chemical potential, fluctuations in the occupations of 
various states are only weakly correlated. Therefore, the fluctuation in the total number of 
particles not in the ground state is microscopic. Hence, a macroscopic condensate fluctuation 
would mean a macroscopic density fluctuation. Even if the particles interact weakly, this 
would mean a macroscopic energy fluctuation. The consequent macroscopic rise in free 
energy would suppress the fluctuation. (See Appendix B for a more formal sketch of this 
argument.) Thus, with interactions producing a finite compressibility, the equivalence of 
the three standard statistical ensembles is assured in the thermodynamic limit, and the 
computationally convenient chemical potential can still be used for isolated, large systems 
@. In the context of Bose liquids, the ideal gas is a theoretical curiosity. Large condensate 
fluctuation is only one of several features for which ignoring interactions gives qualitatively 
incorrect results ||. 

This argument does not address the question of what does happen to condensate fluc- 
tuations of the ideal Bose gas. Furthermore, this is not a totally idle or purely theoretical 
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question. In current experimental work on the trapping and cooling of bosonic atoms, there 
is typically no diffusive particle or thermal energy reservoir While the atoms most 

certainly interact, N ^ oo. Hence, one can ask about the system as a whole rather than 
only describing densities (intensive quantities), which are really just sub- volumes in diffu- 
sive and thermal equilibrium with their (much larger) surroundings. For sub-volumes of 
an infinite system, \i and T give an appropriate description. However, for a finite, isolated 
system taken as a whole, which has a greater impact on the condensate fluctuations, the 
particle interactions or the constraint of fixed total N7 The answer depends on the density 
realized in the particular situation. A practical distinction of a gas from a liquid is that 
the density can be easily varied over many orders of magnitude. In the first successful ex- 
periments [|J, there are noticeable effects of interparticle repulsion; and many of the more 
detailed observations currently underway require a mean field (albeit weak) description of 
the interparticle scattering length to reconcile theory with observations. Nevertheless, it is 
possible to imagine approaching Bose condensation with a box or trap so large and density 
so low that the effects of a given inter-atomic interaction, characterized by a fixed scattering 
length, are negligible, even for density fluctuations of order the equilibrium density. (An 
estimate of the requisite relation of the scattering length, trap parameters and density is 
given in Appendix B.) Even though the Bose-Einstein transition temperature decreases with 
decreasing density, the total energy shift due to a weak fixed-strength inter-particle inter- 
action decreases faster. Also, the actual inter-atomic interactions may not serve to stabilize 
anything. Rather, the gaseous state may itself only be metastable ||. In such situations, 
the equilibrium statistics of the ideal gas are certainly a better starting approximation than 
the equilibrium statistics of the interacting system. 

After a summary of a variety of potentially confusing issues (sec. II), a thoroughly 
elementary analysis of the problem (sec. Ill) suggests that the condensate fractional fluctu- 
ations vanish with increasing N, but all other significant grand canonical predictions have 
vanishing corrections. This is also sufficient to establish the equivalence of using either fixed 
T or fixed total E to characterize the system for large N. The proposed picture provides an 
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explicit prediction (sec. IV) for the condensate fluctuation as well as the values of observ- 
ables, e.g. two-level correlations, that are identically zero with a chemical potential but are 
induced by fixing N. (With a natural normalization, such functions are vanishingly small 
as N — > oo.) The results of a numerical evaluation of the canonical partition function and 
related functions for N from 10 2 to 10 6 (sec. V) confirm these predictions. Some obvious 
conclusions are offered (sec. VI), while comments on details of the numerical work are left to 
Appendix A. Appendix B outlines the simple estimate of the condensate fluctuation damp- 
ing due to repulsive interactions, which allows a comparison with the effect due to fixing 
N. 

II. POTENTIAL ISSUES 

It is only the non-interacting particles in the ground state of a trap or confining potential 
that do not satisfy the hypotheses of the standard demonstration || of the equivalence of the 
grand canonical and canonical ensembles in the thermodynamic limit. Hence, the questions 
raised here only arise if the ground state occupation is macroscopic. At ultra-low T when 
almost all particles are in the ground state, the condensate serves as a particle reservoir for 
all the excited states, and so some form of the grand canonical description for excited states 
should be valid in that domain. But what about intermediate T's? Is the inequivalence 
of chemical potential and fixed N limited to the size of the ground state fluctuations? If 
the condensate manifested the boson propensity for large fluctuations and there were any 
macroscopic fluctuation in the condensate number, it would have to be accompanied by 
correlations between the various occupation numbers. (Such correlations are identically 
zero for the grand canonical ideal gas.) There need not be any macroscopic fluctuation in 
the average density because the total number is fixed. Yet, larger than anticipated exited 
state fluctuations and correlations might lead to larger fluctuations in the total E at fixed 
T. And were this the case, the equivalence of fixing E and fixing T might be lost in the 
thermodynamic limit. 
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Chemical potential is not just a calculational convenience. There is really no practical 
alternative for analytic calculations because not much is known directly about the large but 
fixed N asymptotics of the canonical or microcanonical partition functions, even for systems 
as simple as the ideal Bose gas. If this analytic tool were lost, theory would be reduced 
almost entirely to numerical techniques. 

III. FIXED- N STATISTICS 

The resolution of these conundrums lies in the observation that the grand canonical ex- 
cited state occupations in the thermodynamic limit are independent of not only the conden- 
sate fluctuations but the condensate occupation itself. Hence, if the behavior of the excited 
state occupancies can be reliably estimated using the concept of a chemical potential, one 
can deduce the behavior of the condensate from the constraint of fixed N. This argument is 
really just a minor extension of the traditional one used to compute the condensate fraction 
0|T|]. In particular, it goes as follows. 

Let i label the one-particle (or trap) states and E{ be their energies. Take i = to be the 
lowest energy level, and take Si = 0. In the presence of a chemical potential /i, the mean 
occupation numbers iVj for non-interacting bosons are 

Ni = e ( £i -M)/T _ I ■ W 

With the chosen zero of energy, 

e " M/T = 1 + TT 
N 

- A- 1 (2) 

(defining the fugacity A, to be used later). Once N ^> 1 (which may still be for Nq -C N), 
the explicit fixed-T N dependence of N i>0 is O(l/N ). The expression for the expected 
total number of particles with i > 0, N e , and how it depends on \x is determined by the 
density of states. For an isotropic harmonic oscillator potential in three dimensions with 
level spacing e, 
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N e = C(3) (T/e) 3 



(3) 



as long as N e < N and T/e ^> 1 ||. Under the latter condition, the asymptotic behavior of 
the sum over states is given by an integral. (C(3) ~ 1.202 is the Riemann Zeta function.) 
Under these circumstances, the fixed-T corrections to eq. (0) are 0(1/ Nq). The root-mean- 
square fluctuation of any occupation number is precisely 



so AN e /N e ~ 0(1/ y/Wl). The corrections to eq. (j5|) for \i not exactly zero are again 



The success of using a fi to characterize a system with a large but fixed total number 
of particles N relies on the fact that each individual energy level is a system in diffusive 
equilibrium with the much larger remainder of the total system. This remainder acts as the 
single level's particle reservoir. Once Nq is not much less than N, the utility of /i is no longer 
clear. Certainly there exists no yet-much-larger particle reservoir for the ground state. 

Referring back to eq. ([1]), once N is large, the only role of the particular value of \i is 
to determine N . The A^ >0 are insensitive to /i or No. So, if we consider each individual 
excited level with i > as a system in contact with the reservoir of all the other levels, we 
need not know exactly what the chemical potential actually is, only that it is nearly zero. 
In fact, there need not be any precise meaning to /i, only that it is nearly zero. It may be 
impossible to disentangle the effects of "/i ^ 0" from other 1/N consequences of fixing the 
total N. From this perspective, No is determined not by a /i but by N and N e : 



AN, = ^Ni(Ni + 1) . 



(4) 



For the isotropic oscillator, this implies 




(5) 



O(l/N ). 



N = N-N e . 



(6) 



However, this is precisely the same value of iVo that is deduced from eq. ([I]) when N is 



interpreted as an expectation in the presence of an external //. 
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At the level of occupation expectations, the assignments given by eq. ([!]) for i > 
minimize the total free energy (energy minus Tx entropy) irrespective of the actual value of 
iVo or iV as long as N e is fixed. This is because adding or removing particles from the i — 
condensate changes neither the energy nor the entropy of the entire system. Hence, for large 
N , the occupation numbers for % > are unchanged from their grand canonical values if, 
instead of being determined by a diffusive equilibrium, N is fixed at some value and Nq is 
large. Once there is a condensate, the only thing that can change as particles are added at 
fixed T is N . 

The total expected energy (E) at fixed T depends only on the i > occupations. Thus 
canonical and grand canonical evaluations of the total energy must agree as N — > oo. For 
the isotropic harmonic trap 

(E) = |W 3 (7) 



7T 4 



-TN 



30C(3) 

Since it is a canonical ensemble identity that the root-mean-square total energy fluctuation 
satisfies 



AE=^T^, (8) 

the equivalence of the canonical and microcanonical ensembles is assured as long as N e — > oo 
because AE/E ~ 0{1/ y/W e ). (This is true for any trapping potential, not just the explicit 
example given.) 



IV. FLUCTUATION ESTIMATES 

From the discussion above, it is expected that all occupations approach their grand 
canonical values as — > oo, even if either or both and E are fixed. One can go further 
and estimate the leading behavior of various quantities that vanish in this limit. As examples 
I consider the condensate fluctuations and the occupation correlations between levels. 
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As long as iVo <C N, the root-mean-square fluctuation in the condensate number, ANq, 
satisfies eq. (|). Once N ~ O(N), eq. (j|) implies 

AN = AN e . (9) 

The cross over between these two behaviors is an example of the phenomena that make a 
direct analysis of the fixed- N partition function difficult. It is appropriate to introduce the 
"critical" temperature T c , given by the point at which N e reaches iV or, rather, at which Nq 
goes from macroscopic to microscopic. For the isotropic harmonic potential, eq.@ implies 

T c = N^ 3 C{3)- 1/a e . (10) 

As TV increases, T c remains fixed in absolute, physical units only if the trap size is increased, 
e.g. e decreased. The transition occurs when the central density in the trap reaches the 
infinite volume critical value ||. In terms of the natural temperature variable for the study 
of Bose-Einstein condensation, T/T c , the transition between eq. (f|) and eq. @ takes place 
in a vanishingly small interval as iV — > oo. 

In the thermodynamic limit with N, N , and N e all very large, eqs. (|5||5|,|6|,P,|TU|) can be 
combined to give a simple estimate of the leading behavior: 



AN 1 (T/T c ) 3 / 2 / 7T 2 \ 1/2 
1% ~ Jfil-{T/T c y V6C(3); 



(11) 

A set of quantities that are of interest in the calculation of the angular dependence of light 



scattering off cold, trapped atoms |]T0] are the two-level occupation expectations, (rhiiij). (I 
use the notation "n," for the actual ith level occupation number in a particular configuration 
of the thermal ensemble.) In the grand canonical analysis of an ideal Bose gas, these are 
given precisely by NiNj. In particular, there is no correlation between the fluctuations in 
one level and another. However, with iV fixed, this cannot be exactly true. A refinement 
of the argument of the previous section allows one to estimate the leading behavior of these 
correlations. As an example consider the two states with the largest fluctuations, % = and 
j = 1, because their fixed- iV induced correlation must, therefore, be the largest: 
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At fixed N, if no fluctuates down, say, then n e must fluctuate up by an equal amount. 
The impact on the n i>Q can be estimated by computing the particular expected iVj given 
that N e is larger than its original equilibrium value by the negative of the i — fluctuation. 
This implies (writing Arij for rii — (rii)) 

^HA^HA.l)(-^) t| (12) 

The fugacity, A, is defined by eq.(0). For the isotropic, harmonic trap in the thermodynamic 
limit, this can be evaluated to give (with the natural normalization factor NqNi) 



V. NUMERICAL EVALUATION OF THE CANONICAL ENSEMBLE 



The canonical partition function, Z(N, T), of a trapped, ideal Bose gas can be represented 



as 



Z(N,T) = — dz e iNz Y[ e- nm£m/T ~ inmZ (14) 

71 ^ _7r m=0n m =0 

where n m is the number of particles in the state labeled by m with energy e m . The integral 
over z implements the constraint N = n m . For the isotropic harmonic potential in three 
dimensions, it is convenient to let m label the energy levels e m = me, with the associated 
degeneracy of \{m + l)(m + 2) for m — 0, 1, 2,.... The infinite sums over occupations can be 
done explicitly. Occupation expectations and correlations can be represented similarly by 
simple modifications of the integrand, i.e. extra weight factors of or UiUj. If one truncates 
the infinite product over energy levels m at some finite M max , this yields a form that can be 
evaluated numerically. One can study the convergence in M to test whether the asymptotic 
values of thermal expectations have plausibly been reached. [Useful numerical strategies 
and some details of the evaluations are provided in Appendix A.] 

Fig. (|I|) shows the results of calculations of N . In particular, the solid lines are the 
numerically computed values of N /N for N = 10 2 , 10 3 , 10 4 , and 10 6 , plotted versus T/T c , 
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where T c is given by eq. (|10"D appropriate to each N . The dotted lines are grand canonical 
predictions for iV = 10 2 (small dots) and the iV — > oo limit, 1 — (T/T c ) 3 (large dots). Note 
that the grand canonical predictions were computed as sums over states using eqs. (0,0) and 
involve no approximations regarding N. The comparison of the two statistical ensembles 
is made by identifying the value of the grand canonical (N) with the precise canonical N. 
The canonical numerical calculations clearly approach the N — > oo grand canonical form as 
a limiting value with increasing N. For intermediate values of T/T c , e.g. 0.6, the fractional 
discrepancy between the canonical N and iV — > oo , i.e. difference-divided- by- value, appears 
to be decreasing roughly like jV~ a33 . 

The differences between canonical and grand canonical values for N are displayed in 
another way in fig. (0). The fractional discrepancy between the two evaluations are 
plotted for N = 10 2 , 10 3 , and 10 4 versus T/T c . Here, "fractional discrepancy" means 

^grand canonical _ ^nonicaly N grand canonical _ ^ ^ ^ evaluations giye Nq / N 

very near to 1. So the ratio plotted in fig. (0) plummets, but it is not an effective way 
to characterize the difference between fixed N and fixed /i. (For that region, a more in- 
formative variable would be iV\.) For intermediate values of T/T c , the curves of fig. (0) 
decrease roughly like N~ 1A5 . So, not only does the canonical N approach N(l — (T/T c ) 3 ) 
as N — > oo, it does so approximately as predicted by the simple grand canonical calculation. 
It is the next correction, the difference between the two ensembles' predictions at a given 
iV (as illustrated in fig. (Q)) that reflects the residual difference in physics between the en- 
sembles. This difference is particularly pronounced as A^ makes the transition from micro- 
to macroscopic just below T c . There, the grand canonical - canonical discrepancy decreases 
only very slowly with N. The sign and shape of the difference is such that the canonical iVo 
does not rise quite as sharply as the grand canonical, but the width of the relevant region 
of T/T c vanishes with increasing N. Above T c , the distinction between fixing N and fixing 
H has rapidly vanishing consequences. 

The dashed lines in fig. ([I]) are the results of a numerical evaluation of the canonical 
AiVo/[iVo(iVo + 1)] 1/2 versus the same T/T c 's for N = 10 2 , 10 3 , and 10 4 . For T > T c , this 
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ratio approaches 1, in agreement with the grand canonical eq. (f|). However, for T < T c , it 
goes to zero, more dramatically with increasing N. This same AN data is plotted again 
on a log scale as the solid lines in fig. (|3|). The dotted lines are plots of eq. ([□]) for the 
same N's. As long as neither Nq nor N e are too small, eq. (|TT| ) clearly captures the N and 
T dependence of AN , and the agreement improves with increasing N. In particular, the 
fractional discrepancy between the canonical and eq. ( pTT|) values appears to go roughly like 
N-°- 25 . 

The canonical, normalized, fluctuation correlation, —(AnoAni)/NoNi, is plotted (solid 
lines) on a log scale versus T/T c for N = 10 2 , 10 3 , and 10 4 in fig. (f|). The overall minus sign 
is because the correlation is, indeed, negative. The dotted lines are eq. (p~3|) for the same 
three iV's, and again the agreement improves with N; this time the fractional discrepancy 
appears to go roughly like iV~ a33 . 

The discrepancies between the numerical evaluations and the simple formulae are largest 
for T's such that either iVo or iV e are not very large. These are vanishingly small intervals 
of T/T c for N —> oo. 

The expected i = 1 occupation, Aq, was evaluated to prepare fig. (Q). The agreement 
with eq. ([I]) with /i = was such that the leading fractional discrepancy was accounted for 
by just the leading l/N correction already included in eq. (P, i.e. T/Nq. 

The particular computer code used for the results presented was checked against hand 
calculations for small N. For large N, a criterion for validity was stability under changes in 
the several parameters that should not effect the final answers. Eventually, at high enough 
iV (different values for different observables) the ranges of stability in these parameters 
shrunk to zero. Typically, the practical limitation was the digits of precision available for 
intermediate results. The code was written to evaluate Nq below T c , and specifically for 
Nq plausible results were obtained for much higher N than presented. No effort was made 
to modify the numerical strategy to facilitate calculation of the other quantities discussed; 
presumably those calculations could be extended to higher iV with algorithmic improvements 
that avoided the simultaneous evaluation of numbers of vastly different magnitudes. 
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VI. DISCUSSION AND CONCLUSIONS 



The general arguments presented here, while heuristic, have an internal consistency. For 
example, to compute AN e , which is used implicitly in eqs. ( |12| , |T3|) , one assumes that the 
correlations between level occupations are negligible. One then deduces non-zero correlations 
that are induced by particle conservation. However, the induced correlations are, indeed, 
small enough to be ignored in the calculation of the leading behavior of AiV e and of the 
correlations themselves. 

This is nowhere near to a "theory" of the large N asymptotics of the canonical ideal Bose 
gas. The leading behavior of some interesting observables were estimated and confirmed 
numerically. But in these cases, the leading behavior either was simply given by or could 
be deduced from the grand canonical ensemble. The next level of approximation, e.g. to 
account for fig. (FJ), would require a detailed analysis of the canonical or microcanonical 
partition function and may be very difficult to determine analytically. 

Starting with the grand canonical description with /x and T as independent variables, 
one finds large fluctuations in N below T c . Hence, fixing N may have been expected to be of 
some consequence. However, the grand canonical total energy fluctuations are always small 
and vanish relative to the mean total energy in the thermodynamic limit. Nothing special 
happens in E at T c . So fixing E should have no dramatic consequences. Overall, the switch 
from T to E should be of even less consequence than the switch from /j, to N. A direct 
numerical evaluation of the microcanonical partition function would be considerably more 
difficult. 

However, from a practical standpoint, the modest results here are useful. The largest 
consequence of going from a chemical potential to fixed N for an ideal Bose gas is that 
the ground state number fluctuations are always microscopic; the leading behavior of all 
expected level occupations are unchanged. This is sufficient to further imply that fixing 
the total E is no different from the analytically simpler fixing of T in the thermodynamic 
limit. The leading behaviors two- level expectations, (niUj) for i ^ j, are unchanged because 
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the induced correlations vanish as iV — > oo. For large, fixed N, the corrections to these 
behaviors are unlikely to be of any practical importance. As discussed in Appendix B, for a 
gas with replusive interactions, the consequence of fixing iV dominates over the interaction 
effects in damping the ground state number fluctuations only if the pairwise energy in the 
ground state is less than 0{N~ 2 ^e). 
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APPENDIX A: NUMERICAL STRATEGIES 



For the isotropic harmonic potential in three dimensions and a maximum energy level 
M max , eq. (|i"4l) takes the explicit form 

1 p-n A^max / i \ 2 (m+1) (m+2) 

J ^ 7T m=0 v 7 



A rather primitive C program on a Sun SPARC10 for integrating eq. fllBD and related 
functions was sufficient to generate the numerical results presented in the figures, with the 
size of N limited by the use of double-precision arithmetic. A few general observations may 
prove to be of some value in any future effort to perform comparable calculations. 

Instead of simply truncating the product over energy levels m at some large value M max , 
one can use Maxwell-Boltzmann statistics for all levels m > M max and derive an approximate 
closed form for the contribution to the integrand of all levels above M max . This vastly 
improves the rate of convergence in M max because for modest m's (e.g. 6 x T/e) there are 
still quite a few particles at that m or higher, but the occupations of individual states are 
rarely greater than 1. 
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By far the most rapid variation of the integrand for large N comes from the factor e . 
The integration algorithm should reflect this knowledge. For example, one can divide z into 
intervals of n/AN and integrate each interval accordingly. (For the largest of iV's it proved 
sufficient to take a single point in each such interval.) 

An overall factor in Z has no effect on physical observables. This can be used to consider- 
able advantage. Here are a couple of examples: One can evaluate the products of very large 
numbers logarithmically, i.e. sum the phase and log(modulus) of the various complex fac- 
tors. An overall shift before exponentiation and addition (integration) keeps numbers from 
getting too big. Also, observables are independent of shifts of the whole energy spectrum 
by the ground state energy e . It is convenient to take this non-zero to check the numerical 
independence. Taking Eq ^ can also dramatically alter the character of the integrand of 
eq. (p~5|) — note the (analytically integrable) singularity at z = for e = 0. 

It is, of course, sufficient to integrate only < z < ir. With suitable choice of So, starting 
at z = one can integrate outward, test the convergence, and exit the integration long 
before reaching z = tt. 

APPENDIX B: INTERACTION DAMPING OF GRAND CANONICAL 
OCCUPATION FLUCTUATIONS 

The effect of a weak repulsive interaction on condensate fluctuations can be estimated as 
follows. Let uq represent the number of particles in the ground state. The leading effect of a 
weak, pairwise repulsion at low T, when most of the particles are in the ground state, is to 
raise the energy of those particles from (a convenient n-independent normalization of the 
non-interacting ground state energy) to Auq, where A is the positive two-particle interaction 
contribution to the ground state energy. In natural oscillator units (h — m — — 1), A 
is related to the conventionally defined scattering length a by A = a/^/2ir, at least if the 
interaction effects are weak enough to be treated in mean field theory. Focus on the terms 
in the grand canonical partition function that refer only to the ground state: 
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Qo(fi, T) = J2 e^ T - Xn o/ T . (16) 

n =0 

Unlike the A = case, one can now get large N (= (no)) with \i > 0. Then, the sum can be 
considered as an integral over no, whose integrand is a Gaussian peaked at n = N = /x/2A 
with width AiV = y/T/X. Hence, AN /N = y/T/(XN§), in contrast to the A = situation, 
in which AN /N is 0(1). So grand canonical condensate number fluctuations are small if 
the interaction contribution to the ground state energy is large compared to the temperature. 

If sufficiently strong, interatomic repulsion will certainly be effective at damping conden- 
sate fluctuations at fixed N, giving AN ~ y^TJ\. This effect will dominate (i.e. enforce 
smaller fluctuations) over the non-interacting AN ~ a/ (T/e) 3 estimated in sections III and 
IV when A/e > (T/e)- 2 . (e is the trap level spacing, and T c /e ~ N 1 ^ 3 .) 
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FIGURES 

FIG. 1. Canonical N /N for N = 10 2 , 10 3 , 10 4 , and 10 6 (solid lines), grand canonical N /N for 
N = 10 2 (small dots), the grand canonical N — ► oo limit (large dots), and the normalized canonical 
condensate RMS fluctuations (dashed lines) for N = 10 2 , 10 3 , and 10 4 vesus T/T c . 

FIG. 2. Comparison of the canonical and grand canonical values for Nq as fractional discrep- 
ancies on a log scale for N = 10 2 , 10 3 , and 10 4 versus T/T c . 

FIG. 3. A log plot of the canonical condensate RMS fluctuations (solid lines) and the simple 
eq. (0) estimates (dotted lines) for N = 10 2 , 10 3 , and 10 4 . 

FIG. 4. A log plot of -1 times the normalized, canonical 0-1 level correlations (solid lines) and 
the simple eq. (|l^ ) estimates (dotted lines) for A^ = 10 2 , 10 3 , and 10 4 . 
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fig. 1 -- Condensate fluctuations ... (H.D. Politzer) 




fig. 3 — Condensate fluctuations ... (H.D. Politzer) 
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fig. 4 - Condensate fluctuations ... (H.D. Politzer) 



